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ABSTRACT 

We  study  multivariate  B-splines  M  obtained  as  "shadows"  of 

parallelepipeds,  and  spaces  spanned  by  their  translates  S  =  span  M(  •  -j). 

m 

jSZ 

Recurrence  relations  for  M  are  obtained  and  a  necessary  condition  for  the 
stability  of  the  B-spline  basis  is  given.  We  further  determine  the 
polynomials  contained  in  S  and  the  optimal  degree  of  approximation  from  S. 
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SIGNIFICANCE  AND  EXPLANATION 

Local  support  bases  for  piecewise  polynomial  spaces  are  important  for 
applications  such  as  finite  element  methods,  data  fitting  etc.  In  [BH1 ]  a 
general  construction  principle  for  such  "B-splines*  was  described.  A  special 
case  are  the  sc  called  box-splines.  They  have  a  particularly  regular 
discontinuity  pattern  and  coincide  in  special  cases  with  standard  finite 
elements. 

It  is  hoped  that  using  translates  of  box-splines  will  lead,  at  least  in 
two  variables,  to  a  unified  theory  for  piecewise  polynomial  functions  on 
regular  meshes. 

This  note  is  a  first  attempt  in  this  direction  and  deals  with  basic 
approximation  properties  of  translates  of  one  box-spline  such  as  stability, 
degree  of  approximation  etcv 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


B-SPLINES  PROM  PARALLELEPIPEDS 


C.  de  Boor1  and  K.  Hollig'*^ 

O.  Introduction.  Following  [BH^l,  we  define  the  B— spline  M^  as  the  ■— shadow  of  the 
polyhedral  convex  jody  B  C_  Rn  ,  i.e.,  as  the  distribution  on  R®  given  by  the  rule 

(1)  Mb*  :  =  /B  $°P  .  all  i  S  D(R*)  . 

Here,  P:Rn — *  g":*i — ►  (x(i))™  is  the  canonical  projection  and  /y  denotes  the 
);-dimensional  integral  over  K  in  case  dim  K  =  k  ,  i*e.,  K  spans  a  k-dimensional  flat* 
It  is  obvious  that  MB  is  nonnegative,  with  supp  Mfi  P(B]  •  It  is  easy  to  see 

that  MB  is  a  piecewise  polynomial  function  of  degree  *  n-m  once  one  knows  the  recurrence 
relation  [BH1 ] 

(2)  DpzMB  =  -Ii<z|ni>  Mb  ,  all  7.  6  Rn  . 

i 

Here, 

Uyf  :»  t  y(i )  Dif  , 

with  D^f  the  partial  derivative  of  f  with  respect  to  its  i-th  argument.  Further, 
denotes  the  typical  (n-1 (-dimensional  polyhedron  of  which  the  boundary  of  B  consists,  and 

n^  denotes  the  corresponding  outward  normal.  Finally,  <*|*>  denotes  the  scalar  product. 

In  principle,  KB  can  be  evaluated  with  the  aid  of  the  stable  recurrence  [BH1  ) 

(3)  (n-m)M-(Pz)  -  Z . <b  -z | n, >  M  (Pz)  ,  all  z  6  Rn  , 

°  111 

with  b^  an  arbitrary  point  in  the  flat  spanned  by  B^  . 

Cases  of  particular  interest  are: 

(l)  the  simplex  spline,  obtained  when  B  is  a  simplex.  These  B-splines  were 
introduced  in  tB)  following  up  on  (S)  and  have  already  been  studied  intensively,  mostly  by 

W.  Dahmen  and  C.  A.  Micchelli  lM1_2J,  [D, _4J ,  lDM,_3],  but  alo  by  Goodman  &  Lee  [GL] , 

Hakopian  [Hkj_3),  and  by  H&llig  [Hl^^]  . 
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(ii)  the  truncated  power3  or  cone  splines,  obtained  when  B  is  a  proper  polyhedral 
cone  spanned  by  some  basis  for  Rn  .  These  were  introduced  by  Dahmen  [Djl  and  further 
studied  by  Dahmen  and  Micchelli  in  [DMj]  .  For  example,  they  show  that,  near  an  extreme 
point  of  its  support,  a  simplex  spline  coincides  with  a  truncated  power. 

(iii)  the  box  spline,  obtained  when  B  is  a  parallelepiped.  These  splines  were 
introduced  in  [BD]  and  are  the  object  of  study  of  the  present  note. 

To  be  precise,  a  box  spline  is,  by  a  definition  slightly  more  general  than  the  one  in 

(BD) ,  a  distribution  M_  on  R®  given  by  the  rule 

r 

(4)  M-i*  I— ♦  /  ♦(  £  X(i)5  )  dX 

[0,1)r  i-1 


for  some  sequence  H  :a  (t  .  If  dim  <E>  =  m  ,  then 

=  MD  /  vol  B 

a  B  r 

with 

r  « 

B  :=  {  I  X(i)C  :  X  6  [ 0, 1  ] r) 

i»1  „ 

the  parallelepiped  spanned  by  some  linearly  independent  sequence  in  Rn  for  which 

P^  **  Cj,  ,  all  i  .  We  would  like,  though,  bo  consider  M_  also  in  case  dim  <S>  <  m  . 

For  this,  we  find  it  convenient  to  enlarge  the  above  definition  of  the  B-spline 
by  allowing  P  in  (1)  to  be  an  arbitrary  linear  map  on  B  into  R®  .  Then  (1)  defines 
Mg  as  the  P -shadow  of  B  .  One  checks  that  this  leaves  the  recurrence  relations  (2)  and 
(3)  unchanged  (see  Sect.1). 

In  these  terms,  the  box  spline  M_  defined  by  (4)  is  the  P-shadow  of  the  box  (0,1 ]r, 
with  P  the  linear  map  given  by 

r 

PX  i=  £  X(i)5.  . 
i-1  1 

Here  is  an  outline  of  the  paper.  We  discuss  P-shadows  in  Section  1.  In  Section  2,  we 
give  some  basic  information  about  the  box  spline  M-  ,  such  as  its  recurrence  relations, 
its  Fourier  transform,  and  its  relationship  to  the  difference  operator  d_  and  to  the 
truncated  powers.  We  show  in  Section  3  that  it  is  usually  possible  to  make  a  partition  of 
unity  out  of  the  box  spline  and  certain  of  its  translates  in  many  ways.  We  use  this  fact  in 
Section  4  to  show  that  the  box  spline  and  its  translates  are  usually  globally  linearly 


dependent,  thus  destroying  all  hopes  for  stability  or  the  existence  of  a  set  of  dual  linear 
functionals  for  such  sets  except  in  special  circumstances.  One  such  is  discussed  in  [BH2]. 
In  the  remaining  sections,  we  consider  the  space 

S_  :»  {  Iy  a(j  >N  :  a  @  RV} 

wi  th 

M j  :=  ( • -j )  ,  all  j  6  V  :*  a™  , 

„  „  m 

and  under  the  assumption  that  a  C  V  and  that  <a>  =  R  .In  Section  5,  we  determine  all 

polynomials  in  S=  as  well  as  the  largest  k  for  which  all  polynomials  of  (total) 

degree  k  or  less  are  contained  in  S_  .  We  use  this  information  in  Section  6  to  construct 

a  quasi-interpolant  from  S=  and  thereby  to  obtain  statements  about  the  degree  of 

approximation  obtainable  from  S_  (x  I— *■  f  (x/h)  :  f  6  S_}  as  h  - ►  0  . 

-«h  - 

We  could  have  obtained  our  results  concerning  S_  with  the  aid  of  the  general  theory 
of  spaces  spanned  by  translates  of  a  fixed  function  developped  by  Fix  and  Strang  [FS], 
particularly  if  we  had  been  content  to  discuss  only  1^  .  We  chose  to  derive  our  results 
directly  since  it  seems  no  more  effort  to  do  this  than  it  is  to  verify  that  the  general 
conditions  given  in  (FSl  are  satisfied  for  our  specific  examples. 

We  point  out  in  Section  2  that  S_  £  ,  with 

d  :  =>  max  {  r  :  <a\z>  =  rb  for  all  Z  £  H  with  |z|  «  r  } 

(see  Section  2  for  how  we  treat  the  sequence  a  as  a  set).  This  raises  the  question  of  the 

relationship  of  S_  to  the  space  of  all  pp  functions  in  on  the  same  mesh  and  of 

degree  <  | H | — m  .  We  study  this  difficult  question  in  [BHj]  just  for  m  =  2  and  mainly 
only  for  the  3-direction  mesh,  i.e.,  for  ran  a  »  {e^  e^,  e1+e2^  * 

Rotation.  With  A  £  R®  ,we  denote  by  (A)  the  convex  hull  of  A  and  by  <A>  its 

linear  span.  We  use  x(r)  for  the  r-th  entry  of  the  vector  x  .  For  x  6  R01  and 

j  6  Z+m  ,  the  number  x^  is  computed  as 

xi  x(1)i(,)  ...  x(ra)-* (m) , 

as  usual.  We  denote  by  «  the  class  of  all  polynomials  (on  Rm  ),  and  by  its  subspace 

made  up  of  those  of  total  degree  no  larger  than  k  .  Thu*. 

»  :=  (x  I — ♦  £  a ( j  )x^ } 

*  lil«k 
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and,  more 


with  |j|  :»  j  <  1 )  +  ...  +  j(m)  .  We  also  use  :=  ...  Dm^m* 

generally,  p(D)  !■  E^  a(k)D  in  case  p:x  • — ►  E  a(k)x  .  Here,  we  use  again  the 
notation  D^f  for  the  partial  derivative  with  respect  to  its  i-th  argument  of  the 
function  f  with  domain  in  r“  .  we  also  use  the  notation  :»  E^yti^f  . 

For  a  sequence  of  vectors  in  K°  ,  such  as  H  =■  (C^,...,?)  ,  we  use 


we  also  use  A_  :• 


Finally,  we  denote 


1  r 

A-  ...  A  and  :  *  V  , .  .  V  ,  wi th 
S  4r  a  M  4r 

A  f  s»  f(*+y)  -  f  ,  V  f  :»  f  -  ft*-y)  . 
by  D(K°)  the  space  of  tempered  distributions  on 
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1.  P-shadows.  As  defined  in  the  introduction,  the  P-shadow  of  a  convex  polyhedron  B 
in  Rn  is  the  distribution  M  on  Rm  given  by  the  rule 

M  :  ♦  I — *  /B  *°P  ,  all  ♦  6  D(Rm)  , 
with  P:Rn  — ♦  Rm  an  affine  map. 

We  claim  that  the  recurrence  relations  for  B-spli.^es  established  in  [BH^  ,*  Theorem  21 

remain  valid  for  these  more  general  B-splines  and  state  this  in  the  following  theorem,  for 

the  record.  For  this,  we  make  the  assumption  that  P  is  a  linear  map,  i,e,,  PO  =  0  ,  This 

can  always  be  achieved  by  a  translation  in  Rn  •  Further,  we  assume  that  B  is  proper, 

i.e,,  n-dimensionai.  If  r  :=*  dim  B  <  n  ,  then  this  can  be  achieved  by  restricting  P  to 

the  affine  hull  of  8  and  identifying  this  hull  with  Rr  .  Given  that  B  i3  a  proper 

convex  polyhedron,  its  boundary  is  made  up  of  (n-1 ) -dimensional  convex  polyhedra  ,  with 

corresponding  outward  normals  n^  ,  and  b^  denotes  an  arbitrary  point  in  the  affine  hull 

of  .  Further,  D  stands  for  the  first  order  differential  operator  given  by  the  rule 

r 

Of  :=  £  4.D.f 

i=1  1 

in  case  f  has  its  domain  in  Rr  ,  with 

(*. f)(x)  :=  x(i)f (x)  . 
v  1 

Thus  (Df )  ( x )  =  (Dxf)(x)  and  the  adjoint  of  D  is  -E  • 

Theore,  1.  Let  B  be  a  proper  convex  polyhedron  in  Rn  and  let  M  be  its  P-shadow 
in  under  the  linear  map  P:Rn  — +  Rm  .  Then 

( i )  DpzM  =  -  £^<z|n^>  ,  all  z  G  Rn  . 

(ii )  (n-m)M(Pz)  =  £^<b^-z|n^>  M^(Pz)  ,  all  z  G  Rn  . 

( iil )  DM  =  (n-m)M  -  £^<b^|n^>  . 

Proof*  The  proof  is  a  slight  extension  of  the  arguments  for  Theorem  2  of  [BH1 ) .  The 
proof  of  (i)  only  needs  the  additional  observation  that 

(1)  Dy(*op)  =  (Dpy*)op  . 

This  also  implies  that 

(2)  (t*)(Px)  =  (Dpx<0<Px)  -  (Dx($op))(x)  =  (d($oP)  )  (X)  , 
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of  use  in  a  moment.  The  recurrence  relation  (ii)  follows  from  (i)  and  (ill).  As  to  (ill) 
observe  that 

therefore 

-  (dm)* 

and 

Here,  the  last  integral  in  the 
Thus 

(dm)*  -  (n-m)M*  -  if  D,*.(*op)  , 

i-1  B  1  i 

and  the  argument  now  finishes  as  in  (BHj).  | || 


D  * . 
j  3 


1  +  *.D.  , 
3  3 


/B(j!:iDj*j^OP  =  aM*  +  /B(D*)°P 


/  I  D  *  (*0P)  =  nM*  +  /  D(*op)  . 

B  i-1  1  1  B 

first  line  equals  the  last  integral  in  the  second,  by  (2) 
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2.  Basic  properties  of  the  box  spline.  The  box  spline  M_  defined  in  (0.4)  (with 

r  =  n)  is  a  symmetric  function  of  the  sequence  5  =  (?  )”  .  In  other  words,  M=l  =  M=  for 

any  rearrangement  S'  of  the  sequence  S  .  For  this  reason,  we  find  it  excusable  and,  in 

any  case,  convenient  to  treat  S  in  the  sequel  a3  if  it  were  a  set,  of  cardinality  n  , 

rather  than  a  sequence,  even  though  the  set  (5^:  i=1,...,n}  may  well  have  fewer  than  n 
elements.  Thus,  we  write 

n 

Z  X(C)£  instead  of  £  X ( i ) € . 

CSS  i=1  1 

or 

=  instead  of  (C  1 

v  s(i)y1 

for  the  appropriate  subsequence  s(1),  ....  s(n-1)  of  1,  ...,  n  .  In  the  latter  example, 
this  abuse  of  notation  stresses  the  fact  that  it  doesn't  matter  which  one  of  the  possibly 
several  occurrences  of  the  vector  5  in  the  sequence  5  is  being  omitted. 

It  is  clear  that  M..  is  nonnegative  and  that 

(1)  supp  M_  =  {  £X(C)£  i  X  G  [0,1]”}  . 

Further,  from  (0.4), 

(2)  IM=I  =  1 

as  a  linear  functional  on  CMr”)  ,  Also, 

(3)  M-  G  iff  <S>  =  R®  . 

The  recurrence  relations  of  Theorem  1  for  general  B-splines  simplify  for  the  box 
spline  as  follows. 

Proposition  2.1.  z  »  E_  X(£)£  ,  then 

(4)  DM-  =  Z  X  ( c )  (m-.  P  -  M-v, <•-£))  , 

z  -  =\t 

(5)  (n-m)M-(z)  =«  Z  (  X(5)M-V,(z)  +  (1  -  X  ( £ )  )M-V.  (z-£) )  . 

“  £6=  *  aNfe 

Proof.  The  typical  facet  (i.e.,  (n-1 )-dimensional  face)  of  B  »  (0,1 )n  has  the  form 

:=  {u  G  [0,1]=:  M(t>  =  0} 

or  else  the  form  ,  for  some  £  6  d  •  Further,  and  have  the  outward 


normal  -e^  and  e^ 


respectively.  Thus 


<K|  v 


-M(5)  ,  Bi«  B? 
U(C)  ,  Bi=  e^+B^ 


and  (4)  and  (5)  now  follow  from  Theorem  1.(i)  and  (ii) 


Smoothness.  We  associate  with  E  the  number 
(6)  d  :=  max  {  r  :  <S  \  Z>  =  R™  for  all  Z  C  5  with  |z|  »  r  } 

and  say  for  short  that  E  is  d-spanning.  (We  take  d  =  -1  in  case  <=>  *  Rm  ).  The 
number  d  is  of  interest  here  since  it  follows  from  (4)  and  (3)  that  all  r-th  order 
derivatives  of  M._.  are  in  L—  as  long  as  <H\Z>  =  Rm  for  all  Z  £  S  with  |z|  -  r  . 
Thus 


m_  e 


.  (d) 


,(d-1  ) 


Obviously,  d  cannot  be  bigger  than  j H |  —  m  which  is  the  total  degree  of  the  polynomial 

pieces  of  which  M_  consists.  Precisely,  on  each  connected  component  of  the  complement  of 

{  [£N.z]  +  Z  n  :  H  C  z  ,  <S\z>  *  Rm  }  , 

H 

M-  agrees  with  some  polynomial  of  degree  <  |S|  -  m  . 

■samples,  (i)  For  m  -  1  and  5  *  e,  ,  all  ?  6  =  ,  M=  is  just  the  forward  cardinal 
B-spline,  i.e.,  M_  *  M(* ;0, 1 , . . . ,n)  .  For  m  >  1  and  5  containing  only  e^  effl 

(each  at  least  once),  is  the  tensor  product  of  such  univariate  B-splines. 

(ii)  For  m  -  2  and  |E|  «  n  »  3  ,  with  d  *  1  ,  we  obtain  a  standard  linear  finite 


element. 


(iii)  For  m  »  2  and  =  »  (e^,  e^,  e1+*2»  Cj-e^)  '  is  a  piecewise  quadratic 

function  first  studied  by  Zwart  [Z]  and  independently  derived  by  Powell  (P)  and  Sabin  (PS), 
Its  support  is  shown  in  Figure  1  together  with  its  "mesh",  i.e.,  its  lines  of  transition 
from  one  polynomial  piece  to  a  neighboring  piece.  The  dotted  mesh  lines  occur  in  the  above 
references.  Our  construction  makes  clear  that  they  do  not  appear  in  actuality  since  they  do 
not  lie  in  some  one-dimensional  image  P(FJ  of  some  face  F  of  (0,1 ]n  . 

(iv)  Further  examples  for  S  containing  only  e,,  e2,  e^j  and/or  e^ej  can  be 
found  in  Sablonniere's  study  (SI)  of  smooth  finitely  supported  pp  functions  on  regular 
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Figure 


Support  and  meshiines  for  a  C1  -quadratic  box  spline 


meshes.  The  "generalized  triangular  splines"  of  Frederickson  [F^^J  can  now  recognized 
as  spanned  by  the  box  spline  M-.  ,  with  H  containing  each  of  the  three  vectors  e^,  e ^ * 
and  the  same  number  of  times. 

Associated  difference  operator.  It  follows  from  (4)  that 

DCMS  =  mH\5  -  W'5’  =  Wt  ' 


therefore 

(7) 

In  particular. 


D  M~ 
Z  a 


9  H..  v  „ 
Z  a  \  Z 


for  Z  C 


=  V_5  , 


since  =  6  :=  point  evaluation  at  0  .  Therefore 

(B)  /  M-.D=$  =  [  A„<t»  J  ( 0>  ,  for  all  ♦  e  cl  =  l(Rm)  . 

This  close  association  between  the  box  spline  and  the  forward  difference  operator  brings  to 
mind  the  well  known  association 

/  M(x;tQ,...,tn)  <f<n,(x)/n!  dx  *  [tQ, . . . , 
between  the  univariate  B-spline  and  the  divided  difference. 


The  Pourier  transform  of  M_  is  quite  simple, 


(9) 


fl-(x) 


TT  — — 

•  I  i  £. 


-iC'x 


From  this  we  see  that 


(10>  mh  *  Mz  ■  "huz  • 

SyHetries  and  local  structure.  We  pointed  out  earlier  that  M-  does  not  depend  on 

the  order  of  the  vectors  in  the  sequence  H  .  This  is  due  to  the  fact  that  any  linear  map 

in  Rn  which  permutes  the  unit  vectors  leaves  the  box  [0,1 )n  invariant.  Multiplication 

of  some  of  the  unit  vectors  by  -1  will  change  [0,1 ]n  ,  but  (0,1 ]n  can  be  restored  by  a 

subsequent  shift.  Therefore 
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(11) 


M_(- 


Cgs 


in  case  ag  is  obtained  from  a  by  multiplying  each  C  6  3  by  0(C)  G  {-1,1}  •  A 
symmetry  of  a  different  sort  occurs  when  s'  is  the  image  of  a  under  some  invertible 
linear  map  Q  on  B°  •  Precisely, 


(12)  H_  »  |det  Q|  M^oq  . 

Thi3  implies  certain  symmetries  for  M=  in  case  a  »  Qs  . 

The  box  spline  is  particularly  simple  near  an  extreme  point  of  its  support. 


Proposition  2.2.  _I£  e  is  an  extreme  point  of  supp  M=  ,  then,  in  a  neighborhood 

of  e  ,  M=  agrees  with  some  truncated  power  of  degree  |5|  -  m  .  In  particular, 

M_  ( •  +  e )  is  homogeneous  of  degree  | a |  -  m  near  0  . 


Proof.  We  pointed  out  earlier  that 

supp  M_  ■  {  S  A,C  :  Xe[0,1 ]"}  . 

CGS  4 

Thus  any  extreme  point  e  of  the  (closed)  support  is  necessarily  of  the  form 

e  *  EC 

tez 

for  some  z  C  3  for  which,  further,  there  exists  h  6  R  so  that  <n|C>  >  0  for  all 


C  G  a  with 
0 


0(C)  !• 


-1  ,  Cgz 
1  ,  CgS\z 


Therefore,  from  (11), 

~a 

showing  that  it  is  sufficient  to  consider  the  case  that  e  *  0  and,  for  some  h  G  Rm  and 
all  C  G  3  ,  <n|C>  >  0  . 

In  this  case, 


e  :»  dist(0,  (31)  >  0  . 

Therefore,  for  all  test  functions  ♦  with  supp  ♦  £  B£(0)  !“  hall  of  radius  c  and 
center  0  , 

M  -  /  ♦(£  X(C)C)  dA  -  /  ♦(£  A(OC)  dA  .||| 

[0,1}  R“ 

-10- 


3.  Partition  of  unity.  In  this  section,  we  show  that  appropriate  translates  of  an 


appropriately  scaled  version  of  the  box  spline 

M  s-  M= 

form  a  partition  of  unity.  By  a  standard  argument,  this  implies  that  the  space  spanned  by 
these  translates  can,  at  least,  approximate  continuous  functions  (as  the  mesh  size  is 
reduced  by  scaling). 

Proposition  3.  Suppose  H  contains  the  basis  Z  (for  R0  ).  Then 

(1)  E  m(*  -  £  j(CK)  =  1/|det  Z|  . 

jez  cez 

Proof.  Since  R®  is  the  essentially  disjoint  union  of  the  sets 

(  £  (A  (C )-j (C ))c  r  A  e  [0,1 1Z)  ,  j  e  ZZ  , 
cez 

we  find  that 

£.«(•-  £  jUK)*  -  /  (  /  ♦(E  A(C)C  +  £_  11(C)?)  dA)  dl*  . 

jgxZ  Z  [0,1]  R  Z  -vz 

The  change  of  variables  A  (— »  £  A(C)C  carries  this  to 

l* 

I  n_m  (/„♦(*  +  £=sZy(OC  )  dx/|det  Z\  )  dp 
(0,1]  R 

and  this  equals 

!n  *  voln-m10'1 )nm/ldet  Zl  ‘III 

R 

Corollary.  If  =  C  z"  and  <2>  -  r"  ,  then  £  M(’-j)  ■  1  . 

~  ~  je^" 

Proof.  Let  Z  C  2  be  a  basis  (for  R®  ).  Then 

*  {£,  j(C)C  8  j  6  Zz) 

xs  a  subgroup  of  Z®  and  its  factor  group  G  s=  Z™/A  has  (finite)  order  |det  z|  . 
Therefore 

£  M ( *  — j )  =  E  E  M( • -g-j )  =  | G | / |det  z|  ■  1  .  ||| 

jez"  g6G  jeA 
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4.  Linear  independence  of  translates.  For  any  particular  subset  V  of  R®  ,  we 
consider  the  collection  of  translates  Mv  :=  M(*-v)  ,  v  6  V  ,  of  the  box  spline 

M  M_  . 

Such  a  collection  is  always  (algebraically)  linearly  independent:  Indeed,  if  f 

:  =  £y  a(v)Mv  with  VI  i«  supp  a  a  finite  nonempty  set,  then 

(supp  M  )\U  (supp  M  )  f  0 
v  v«W\w  v 

for  some  w  6  w  ,  hence  f  j*  0  . 

We  are  interested  in  considering  nontrivial  sums  of  infinitely  many  translates.  Por 
this,  we  make  the  assumption  that  V  has  no  finite  limit  points.  Then  only  finitely  many 
of  the  translates  have  any  particular  point  in  their  support  and  thus,  for  a  6  Rv  with 
suitably  controlled  growth  at  infinity, 

f  Zy  a(v)  Mv 

defines  a  distribution  on  R®  . 

Assume  that  M  is  a  function,  i.e.,  that  H  contains  a  basis  (for  R®  ).  Then 

SH,v  span  (Mv^V  {ZV  a(v)  Mv  1  a  8  *V} 

is  a  space  of  piecewise  polynomial  functions,  possibly  quite  smooth,  and  it  becomes  of 

interest  to  find  out  to  what  an  extent  (Hyjy  is  a  basis  for  this  space  or  one  of  its 
subspaces.  We  call  (Mv)v  (globally)  linearly  independent  if  the  linear  map 
( 1 )  a  I — *•  Eva(v)  Mv 

is  1-1  on  RV  .  Such  linear  independence  is  a  first  necessary  condition  for  other 

properties  of  interest  to  hold.  One  such  property  is  stability  of  the  basis  (My)v  for 

S_  ,  i.e.,  the  property  that  the  map  (1)  is  bounded  and  bounded  below  on  lm (V)  into 
a  ,V 

I^(Rm)  .  Another  is  the  possibility  of  interpolation  from  S_  y  ,  i.e.,  the  existence  of 
points  pv  ,  typically  with  Mv(py)  j*  0  ,  all  v  6  V  ,  so  that,  for  any  function  f  in 
some  class  K  ,  there  exists  one  and  only  one  s  6  S_  ,.nK  which  agrees  with  f  at 
( Pv > v  *  We  malce  clear  below  that  this  (global)  linear  independence  is  usually  not  present. 
Yet,  as  is  pointed  out  in  (BD),  _if_  Zv  a(v)Mv  »  0  and  a  j*  0  ,  then  there  exists  r  >  0 
so  that,  for  all  v  6  V  ,  a  changes  sign  on  V  '|Br(v)  .  This  implies  that  the  map  (1)  is 
1  - 1  on  *  |  v  . 
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The  space  S_  y  =  span  (MVJV  becomes  interesting  when  V  is  related  to  S  .  By 
assumption,  M  is  a  function,  i.e.,  H  contains  a  basis  Z  (for  R®  ).  Therefore, 
according  to  Proposition  3.1,  the  collection  Mv  ,  v  6  Vz  :=  {  j(CK  :  j  S  ZZ}  forms  a 
partition  of  the  constant  l/|det  z|  .  This  suggests  consideration  of  V  of  the  form  Vz 
for  some  basis  Z  in  H  .  We  go  one  step  further,  though,  and  consider  from  now  on  only 
the  following  normalized  situation: 

(2)  a  £  V  -  Zm  . 

This  is  the  same,  up  to  an  affine  transformation,  as  the  assumption  that  V2  —  v  ^or 
bases  Z  in  a  .  We  abbreviate 

S_  :=  S_  m  . 

- 

Proposition  4.  Under  the  assumption  (2),  (M j )y  is  linearly  dependent  unless 
|det  z|  =  1  for  all  bases  Z  C  3  . 


Proof.  By  assumption,  a  contains  a  basis  Z  for  R™  ,  therefore  (|det  Z|m^) 

provides  a  partition  of  unity  as  does  (M^ )y  ,  by  Proposition  3  and  its  corollary. 

now  |det  z|  j*  1  for  some  basis  Z  in  S  ,  then  V  ,  yet 

Ey  |det  Z|  M.  =  1  =  £y  M  .  HI 

Z 


'Z 

If 


Reaark.  It  would  be  nice  to  know  whether  the  converse  of  this  proposition  holds. 
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5.  The  polynomials  In  S-  .  In  this  section,  we  deternine  S-n*  .  This  information 
is  important  in  the  discussion  of  the  degree  of  approximation  to  smooth  functions 
attainable  from  S_  .  We  continue  to  use  the  abbreviations  and  assumptions  introduced  in 
Section  4. 

Lemma  5.1.  *  n  her  0_  -  inter  A_  . 

Proof.  Recall  from  (2.8)  that 

(i_f)(x)  »  /  M  D_f(*+x)  . 

Therefore  *  n her  A_  D  » nicer  n_.  .  For  the  converse,  observe  that  d_f  »  o  implies 
Jm(*-x)  D_f  «  0  for  all  x  .  Since  M  is  nonnegative  and  of  compact  support,  this 

cannot  hold  for  a  polynomial  f  unless  the  polynomial  D_f  vanishes  identically.  | | | 

We  also  note  that  (2.7)  together  with  summation  by  parts  gives 

<1»  Dz(r  a<j)  Wj)  -  I  (*za)(j)  Msvz(*"j)  '  a11  Z  -  3  • 

Theorem  5.  Let  K  f”)  *  her  D  with  :*  :■  Izc  ■  ,  <3\Z>  *  Rm}  .  Then 

Z«3  1  -  ~  - 

(2)  ms.  -  ms  , 

s 


Proof.  Let 

Z  a(j)  «!  p  e  s=  . 

If  Z  £  =  ,  then,  by  (1),  the  polynomial  Dzp  can  be  written 

V"  1  "aW-11  * 

It  also  <HXz>  *  r"  ,  then  supp  M-  has  zero  measure,  hence  the  polynomial  Dzp  must 
vanish  identically. 

For  the  converse  statement,  ws  prove  by  induction  on  k  that 
(3)  S-  3  i.ni  , 

it  being  trivially  true  for  k  «  -1  .  For  the  induction  step,  we  show  now  that 
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(4) 


(4)  p  6  *)tr,K  implies  q  :»  p  -  £  p( j  )M^  0  *){_1r,|C  • 

Since  p  belongs  to  K  ,  so  does  £  p(j)Mj  ,  by  Lemma  5.1  and  (1)  (making  use  of  the  fact 
that  ker  A_  =»  ker  V=  ).  Thus  it  remains  to  show  that  q  0  1  .  This  is  established  once 

we  show  that,  for  any  y  G  Rm  , 

(Dy)kq  -  0  . 

For  this,  we  note  that,  whenever  Z  S  ,  then 

(D  )SD„  =  (d  )8'’  £  a(C)  D_v, 

y  z  y  5e=\z  zse 

(with  y  =  Z ,  a(5)C  ).  Repeated  application  of  this  formula  justifies  the  claim  that 


(d  f 

V  yl 


ZGS  , |z|<k 


i(Z)(o  ) 


k-jz| 


D  +  Z  #  a(Z)D 

2  Z0S,Z0S  ,  |z|»k 


It  follows  that 


(OJN 


Z@S,Z0S  , |z|-k 


a(Z>(o7P  -  Z(’7p)(j)M3XzC-j  >) 


and  this  is  zero  since,  for  each  Z  C_  3  with  Z  £  3  ,  we  have  £  M~^z(*-j)  *  '  ^y  the 
corollary  to  Proposition  3,  while  |z|  =  k  implies  that  V^p  *  Dzp  i3  some  constant, 
since  p  6  'k  , 

It  is  now  easy  to  complete  the  induction  step.  If  p  9  »)[r|K  *  then,  by  (4),  p  G 
S_  +  1°*  '  hence  P  e  s-  bY  induction  hypothesis.  ||| 

Corollary  1.  For  each  k  ,  the  map  T:  p  i — *  Z  p(j)  fT  carries  1-1  onto 

itself. 

Proof.  We  mentioned  already  in  Section  4  that  T  is  1-1  on  *  .  Thus  it  is 


*  C  S_  iff  n  C  K  :■  Iter  D_  • 

JC  —  a  K  —  it 

ZQz 

Further,  the  differential  operator  Dz  decreases  the  total  degree  of  any  polynomial  by  at 
least  | Z |  and,  for  some  polynomials,  by  exactly  |zj  .  Finally, 

d+1  «  min,  \z\  , 

Z8S 

»  _*  .  . 

This  shows  that  c  ker  Dz  for  all  Z  e  s  .It  also  shows  that,  for  some  z  6  =  ,  |Z| 

■  d+1  ,  hence  DZP  ?  0  for  some  p  6  *d+1  •  III 


6.  Degree  of  approximation  from  S_ 


In  this  last  section,  we  discuss  the  degree  of 


approximation  from  S-  to  a  sufficiently  smooth  function,  as  h  — *  0  ■  Here,  h 

-I*1 

indicates  a  scaling  of  the  mesh,  i.e., 

S_  .  :=  {x  I— *•  f (x/h)  :  f  6  S-}  . 

a,h  = 

theorem  6.  lf_  k  <  d  (wi th  d  given  by  (2.6)),  then  there  exists  a  linear 

functional  A  on  x,  so  that  p  =  Z_,  Ap(*+-j)  m.  for  all  P  6  it,  . 

k  13  “ x 

Proof.  By  Corollary  2  of  Theorem  5,  x^  C  S_  ,  while,  by  Corollary  1  of  Theorem  5, 

the  map  Tip  i — *  E  p(j)  is  1-1  onto  x^  .  Thus 

p  =  Z  (t  ’p) ( j )  M.  ,  all  p  6  x 
3  3  * 

with  T-1  :=  (t,  )  .  Further,  with  T^  the  shift  by  the  vector  i  ,  i.e., 

'  k 

(t^p] Cx)  :=  p(x+i)  ,  all  x  , 

we  have 

E  E  p(  j+i  )M^  =  Z  p(j)*T(»+i)  =  T^Z  p(j)Mj)  , 

showing  that  T  commutes  with  T^  ,  hence  so  does  T-t  .  This  proves  the  theorem,  with 
Ap  :=  (T-,p)(0)  .  HI 


The  theorem  implies  statements  about  degree  of  approximation  to  smooth  functions  from 
S_  in  the  now  standard  quasi-interpolant  fashion:  Let  K  be  the  class  of  functions  which 
belong  locally  to  some  function  space  Kg  ,  e.g.,  to  L1  .  Extend  the  linear  functional  A 
of  the  theorem  to  a  continuous  linear  functional  W  on  Kq  and  with  support  in  supp  M  . 
The  quasi-interpolant  scheme 

Q:K  - *  S_:  f  I — ♦  Z  uf(»+j)M.. 

then  reproduces  x^  and  is  local.  This  implies  that 

f  -  Qf  =  (f-p)  -  Q(f-p)  ,  for  all  p  6  x^ 

and  that 

|(Qg)(x)|  =  |Epg(*+j  )fMx)  |  «  I  pi  max{  Ig  |  suppM  *  •  Mj(x>  * 

Therefore 


with 


|(f  -  Qf )  <x)  |  <  distp(f,  x^) 


p(g)  |g(x)|  +  III!  max{lg.  ..  I  :  M.(x)  *  0>  . 

I  supp^  J 

This  shows  that  Qf  approximates  to  f  locally  as  well  as  local  polynomial  approximation. 
Here  is  a  particular  result  along  these  lines. 


Corollary.  Let  u  be  an  extension  of  X  to  a  continuous  linear  functional  on 

I^tsupp  M)  and  let  Qj,  :=  SjjQS^jj  with  Q:f  I— ■*  Z  yf(*+j)M^  and  (shf)(x)  :»  f(xh)  . 

Then,  for  f  6  l£1c+1>,  If  -  Qhfl„  =  0(hk+1 )  . 

Sharpness.  The  order  0(h^+1 )  is,  in  general,  best  possible  for  the  approximation 

from  S?  to  smooth  functions#  To  see  this,  choose  Z  6  s  (cf.  Theorem  5)  with  |Z|  ■ 

d+1  and  a  polynomial  p  6  with  D^p  *  1  •  If,  for  some  approximating  sequence  8^ 

fl  S-  .  ,  we  have 
a,h 

a  ■  ..  d+1  . 

,8h  '  plL1CO,HB  ■  o(h  >  ' 

it  follows  from  the  standard  Markov  inequality  for  piecewise  polynomials  that 

(1>  ,DZ'8h  '  VP’L^O,!]"  “  °{h> 

for  any  z‘  C  z  with  |z'|  ■  d  .  Set  z  »  Z\Z'  .  Since  the  support  of  is 

contained  in  hyperplanes,  <DZi3h*  ia  Piecewi*®  constant  on  lines  in  the  direction  z  .  On 

the  other  hand,  Dz,p  restricted  to  these  lines  is  of  the  form  (Dz,p)(x  +  tz)  » 

(Oz,p)(x)  +  t  which  contradicts  O).  ||| 


i 


I 

1 


References 


Ibl  C.  de  Boor,  Splines  as  linear  combinations  of  B-splines,  in  Approximation  Theory 
11 ,  G.G.  Lorentz,  C.K.  Chai  and  L.L.  Schumaker  eds..  Academic  Press,  1976,  1-47. 

[BD]  C.  de  Boor  and  R.  DeVore,  Approximation  by  smooth  multivariate  splines,  MRC  TSR 
#2319  (1981) 

I BH^ )  C.  de  Boor  and  K.  HOllig,  Recurrence  relations  for  multivariate  B-splines,  MRC  TSR 
#2215  (1981);  Proc.  Araer.Math.Soc. ,  to  appear. 

(BH2)  C.  de  Boor  and  K.  HOllig,  Bivariate  box  splines  and  smooth  pp  functions  on  regular 
meshes,  ms. 

(D,)  W.  Dahmen,  On  multivariate  B-splines,  SIAM  J.Numer.Anal.  J_7_  (1980)  179-191. 

[D21  W.  Dahmen,  Multivariate  B-splines  -  recurrence  relations  and  linear  combinations 

of  truncated  powers,  in  Multivariate  Approximation  Theory,  W.  Schempp  and  K.  Zeller 
eds.,  BirkhSuser,  Basel,  1979,  64-82. 

(Dj)  w.  Dahmen,  Konstruktion  mehrdimensionaler  B-splines  und  ihre  Anwendungen  auf 

Approximationsprobleme,  in  Numerische  Methoden  der  Approximationstheorie,  Bd.5,  b. 
Collatz,  G.  Meinardus  and  H.  Werner  eds.,  BirkhAuser,  Basel,  1980,  84-110. 

[D4]  w.  Dahmen,  Approximation  by  smooth  multivariate  splines  on  non-uniform  grids,  in 
Quantitative  Approximation,  R.  DeVore  and  K.  Scherer  eds..  Academic  Press,  1980, 
99-114. 

(DM,)  w.  Dahmen  and  C.A.  Micchelli,  On  limits  of  multivariate  B-splines,  MRC  TSR  #2114, 
1980;  J.d 'Anal. Math.  j!9  (1981)  256-278. 

(DM2)  w.  Dahmen  and  C.A.  Micchelli,  On  the  linear  independence  of  multivariate  B-splines. 

I.  Triangulation  of  simploids,  SIAM  J.Numer.Anal.,  to  appear. 

(DMj)  W.  Dahmen  and  C.A.  Micchelli,  On  the  linear  independence  of  multivariate  B-splines. 

II.  Complete  configurations,  Math, Comp.,  to  appear. 

[F,]  P.  0.  Frederickson,  Generalized  triangular  splines,  Mathematics  Report  7-71, 
bakehead  University,  1971. 

[F21  P.  O.  Frederickson,  Quasi-interpolation,  extrapolation,  and  approximation  on  the 
plane,  in  Proc. Manitoba  Conference  on  Numerical  Mathematics,  Winnipeg,  1971,  169- 
176. 

[FS]  G.  Fix  and  G.  Strang,  Fourier  analysis  of  the  finite  element  method  in  Ritz- 
Galerkui  theory.  Studies  Appl.Math.  _48  (1969)  265-273. 

(Gl>)  T.N.T.  Goodman  and  S.b.  bee,  Spline  approximation  operators  of  Bernstein-Schoenberg 
type  in  one  and  two  variables,  J, Approximation  Theory,  to  appear. 

(Hk,)  H.  Hakopran,  On  multivariate  B-splines,  J. Approximation  Theory,  to  appear. 

(Hkj)  H.  Hakopian,  bes  differences  divisees  de  plusieurs  variables  et  les  interpolations 

multidimensionel les  de  types  lagrangien  et  hermitien,  Comptes  Rendus  Acad. Sc. Paris 
292  ( 1981 )  xxx. 


{Hkj)  H.  Hakopian,  Multivariate  divided  differences  and  multivariate  interpolation  of 
Lagrange  and  Herraite  type,  SIAM  J.Numer.Anal.,  to  appear. 

K.  HOllig,  A  remark  on  multivariate  B-splines,  J. Approximation  Theory,  to  appear. 

K.  HOllig,  Multivariate  Splines,  MRC  TSR  #2188  (1981);  SIAM  J.Numer.Anal.,  to 
appear. 

If 

C. A.  Micchelli,  A  constructive  approach  to  Kergin  interpolation  in  R  %  multivariate 
B-splinea  and  Lagrange  interpolation,  MRC  TSR  #1895  (1978);  Rocky  Mountain  J.Math. 
10  (1980)  485-497. 

[M2J  C.A.  Micchelli,  On  a  numerically  efficient  method  for  computing  multivariate 

B-splines,  in  Multivariate  Approximation  Theory,  W.  Schempp  and  K.  Zeller  eds., 
BirkhSuser,  Basel,  1979,  211-248. 

(Pi  M.J.D.  Powell,  Piecewise  quadratic  surface  fitting  for  contour  plotting,  in 

Software  for  Numerical  Mathematics,  D.J.  Evans  ed..  Academic  Press,  1974,  253-271. 

[PS]  M.J.D.  Powell  and  M.A.  Sabin,  Piecewise  quadratic  approximations  on  triangles,  ACM 
Trans. Math. Software  _3_  (1977)  316-325. 

(Sb)  P.  Sablonniere,  unpublished  notes  1980;  part  of  the  material  has  appeared  in 

reports,  e.g.,  in  De  l'existence  de  spline  a  support  borne  sur  une  triangulation 
equilaterale  du  plan.  Publication  ANO-39,  U.E.R.  d'I.E.E.A.  -  Informatique, 
Universite  de  Lille  I,  1981. 

(Si  I.J.  Schoenberg,  Letter  to  Philip  J.  Davis  dated  May  31,  1965. 

(Zl  P.  Zwart,  Multivariate  splines  with  nondegenerate  partitions,  SIAM  J.Numer.Anal.  _10_ 
(1973),  665-673. 


(HI,] 

(hi2i 

(M,l 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWi •«  Of  Enl.r.Q 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  OOVT  ACCESSION  NO. 

1.  RECIPIENT'S  CATALOG  NUMBER 

4.  title  (end  Subtitle) 

B-SPLINES  FROM  PARALLELEPIPEDS 

5.  TYPE  OF  REPORT  8  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 

S.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR^ 

C.  de  Boor  and  K.  H<^.lig 

8.  CONTRACT  OR  GRANT  NUM8ERf*; 

DAAG29-80-C-00  41 
MCS-7927062,  Mod.  1 

»■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 

610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

fork  Unit  Number  3  -  Numerical 
Analysis  &  Computer  Science 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

See  Item  18  below. 

12.  REPORT  DATE 

February  1982 

IS.  NUMBER  OF  PAGES 

20 

1«.  MONITORING  1GENCY  NAME  A  AODRESSfJf  dll  to  font  (ram  Controlling  Otllco) 

IS.  SECURITY  CLASS,  (ol  thlo  report; 

UNCLASSIFIED 

IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  ( o 1  thlo  Roport) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tho  obttroct  ontmrod  In  Block  30,  It  dllloront  (ram  Report) 

18.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office  National  Science  Foundation 

P.  0.  Box  12211  Washington,  D.  C.  20550 

Research  Triangle  Park 

North  Carolina  27709 

19.  KEY  WORDS  (Continue  on  reweree  aid*  1/  necemeery  end  Identity  by  block  number) 

B-splines 
multivariate 
spline  functions 
degree  of  approximation 

20.  ABSTRACT  (Continue  on  reeeree  elde  II  nacaatary  and  Identity  by  block  number) 

We  study  multivariate  B-splines  M  obtained  as  "shadows"  of 
parallelepipeds,  and  spaces  spanned  by  their  translates  S  =  span  M(-  -j). 

jezm 

Recurrence  relations  for  M  are  obtained  and  a  necessary  condition  for  the 
stability  of  the  B-spline  basis  is  given.  We  further  determine  the 
polynomials  contained  in  S  and  the  optimal  degree  of  approximation  from  S. 

00 


FORM 
I  JAN  71 


1473 


EDITION  OF  I  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh.n  Df  BnlorodJ 


